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SUMMARY 


This report collects, documents, and models data relating the 
expected accuracies of tracking variables to be obtained from the FAA's 
Mode S Secondary Surveillance Radar system. The data include measured 
range and azimuth to the tracked aircraft plus the encoded altitude 
transmitted via the Mode S data link. 

A brief summary is made of the Mode S system status and its 
potential applications for aircraft safety improvement including accident 
analysis. FAA flight test results are presented demonstrating Mode S 
range and azimuth accuracy and error characteristics and comparing Mode S 
to the current ATCRBS radar tracking system. Data are also presented that 
describe the expected accuracy and error characteristics of encoded 
altitude. These data are used to formulate mathematical error models of 
the Mode S variables and encoded altitude. A brief analytical assessment 
is made of the real-time tracking accuracy available from using Mode S 
and how it could be improved with down-linked velocity. 
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INTRODUCTION 

The purpose of this report is to collect, model, and document data 
relating to expected accuracies of tracking variables to be obtained 
from the FAA's new Mode S Secondary Surveillance Radar (SSR) system. As 
compared to the current ATC Radar Beacon System (ATCRBS) , the new Mode S 
SSR provides more accurate range and azimuth measurements of the tracked 
aircraft. Mode S also includes the ability to data link additional 
information between the aircraft and the ground. A procurement is about 
to be released by the FAA to obtain and install over 100 of these ground- 
based tracking systems. It is expected that by 1990, all airspace above 
12,500 ft plus relevant terminal areas will have access to Mode S service. 

The Mode S system has been tested by the FAA at their Technical 
Center, and the test data are summarized in the subsequent chapters. The 
data include measured range and azimuth to the tracked aircraft plus the 
encoded aircraft altitude transmitted via the Mode S data link. The 
accuracy of this data is of importance to various aircraft safety 
improvements, including accident investigation. 

Chapter II gives a brief summary of the status of Mode S and the 
potential applications that it has for improved aircraft safety, including 
accident analysis. This includes a summary of uses of the Mode S data 
link capability. 

Chapter III presents FAA Technical Center flight test results that (a) 
demonstrate the range and azimuth measurement accuracy available from Mode S, 
and (b) describe the characteristics of the measurement errors. This in- 
cludes a comparison to the current ATCRBS radar tracking system used for 
air traffic control. 
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Chapter IV presents a summary of error characteristics of the encoding 
altimeter. This information is necessary to obtain the vertical dimension 
of the recreated aircraft flight profile. 

Mathematical error models of the Mode S system and the encoded altitude 
are presented in Chapter V. Then, Qiapter VI summarizes these results and 
makes recommendations for further research in support of aircraft safety. 

Appendix A presents a brief analytical assessment of the real-time 
tracking accuracy that is available from using Mode S and how it could be 
further improved with down-linked velocity information. 
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II 


MODE S APPLICATIONS 


The FAA published a 20-year plan this past year [1] which describes 
the development and upgrading of many systems affecting air traffic control 
(ATC) . Table 1 is the projected summary of naticial airspace system 
activity through 2000. As can be seen, air carrier, commuter and general 
aviation (GA) aircraft numbers are expected to grow by 42%, 175% and 94% 
in the twenty-year span. Air carriers and commuters will all be transponder 
equipped with probable altitude reporting (Mode C or Mode S) cabability. 

The current 200,000 GA aircraft are about 60% transponder equipped [2], 
and this number is ever increasing. 

One major development that enhances aircraft safety is the design and 
testing of a new ATC surveillance radar system known as Mode S Secondary 
Surveillance Radar (SSR) . This development includes upgrading and replace- 
ment of the ground-based system (re: erred to as the "sensor") and replace- 
ment of the airborne component (referred to as the "transponder"). Mode S 
allows for evolutionary growth ir. the ability to communicate between the 
aircraft and the ground; it .ll' ws the sensor to send and receive unique 
messages from individual ai ft transponders (by using ind''vidual discrete 
addresses) . 

Reference 1 states that "Mode S and data link covera 2 will be provided 
by 1990 above an altitude of 12,500 ft above mean sea level and down to the 
surface of qualifying a rports." Further it states that "to obtain air 
traffic control cleara*'.ce, aircraft will have to be equipped with Mode S by 
1990." Also, the datwi link coverage will be extended from 12,500 ft altitude 
down to 6,000 ft abr /e mean sea level by year 2000. This coverage will be 
provided by 137 Moc-j S sensors by 1990 with an additional 60 by 2000. The 
implication of th- above is that all aircraft, to fly IFR, will be required 
to have Mode S transponders. 
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Table 1. l->tal National Airspace System Activity [1] 



1980 

1985 

1990 

2000 

PERCENT GROWTH 
1C30-2000 

NASP Airports 

3163 

3637 

3631 

4000 

26.5 

Aircraft Operations 
(mi 1 1 i ons ) 

134.1 

179.7 

212.7 

290.0 

116.2 

Itinerant Operations 
(mi llions ) 

74.0 

98.8 

117.3 

159.5 

115.5 

Instrument Operations 
(mil lions) 

38.2 

48.1 

54.2 

65.6 

71.7 

IFR Aircraft Handled 
(mi llions) 

30.1 

37.2 

42.2 

53.7 

78.7 

FSS Operations 
(mi 1 1 i ons ) 

64.3 

83.9 

98.6 

139.6 

117.1 

Dorrtestic Enplanements 
(mi llions) 






Ai r Carrier 

278.1 

380.8 

454.0 

589.8 

112.1 

Commuter 

13.1 

21.8 

30.6 

42.0 

220.6 

Ai rcraft Fleet 
(thousands ) 






Air Carrier 

2.4 

2.7 

2.9 

3.4 

42 

Commuter 

1.6 

2.4 

3.2 

4.5 

175 

General Aviation 

210.3 

254.5 

298.1 

408.5 

94 

Pilots 

(thousands) 






Instrument Rated 

247.1 

309.5 

369.6 

481.9 

95.0 

Total Pilots 

814.7 

891.0 

1050.6 

1331.3 

63.4 


The applications of the Mode S data link are still being developed, 
but the following are considered as strong candidates 13] : 

1. Pilot retrieval of "weather products". (Digital weather maps) 

a. Terminal weather conditions; 

b. Winds aloft; 

c. Hazardous wer'^her conditions; and 

d. Route oriented weather. 

2. Air Traffic Advisory Service (ATAS) - This will give to the 
GA pilot surrounding traffic information and provide backup 
protection to the primary ATC control in dense terminal areas; 

3. Minimum safe altitude warning (MSAW) alerts; 

4. Takeoff clearance confirmation; 

5. Altitude assignment confirmation; 

6. Airport surface clearance; 

7. Aircraft position relative to a sensor. This could be 
considered as an alternate navigation aid; and 

8. Automatic sector handoff. 

In additio , the Mode S transponder is an essential element of the 
TCAS II collision avoidance system now being developed. A long-term 
development also being researched is the cockpit display of traffic 
information (CDTI) which could use either the ATAS service or TCAS system 
as the source of adjacent aircraft state information [4,5]. 

Another potential application » of major concern in this report, is 
improved accident investigation capability. Previous studies [6,7] have 
established that ATC tracking radar data can be used to reconstruct the 
trajectory of a maneuvering aircraft. By taking the time history of the 
aircraft position data as measured by the tracker, smoothing and differenti- 
ating this data, and then by using known kinematic relationships governing 
aircraft flight, a reasonably accurate reconstruction of the profile csn 
be generated. This trajectory reconstruction is valuable for analyzing and 
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deducing the cause of aircraft accidents. This process is limited, however, 
by the accuracy of the tracking data which consists of measurements of slant 
range and azimuth from the tracker s'^te plus encoded altitude. The Mode S 
SSR has two features which enhance significantly the ability to recon- 
struct the aircraft trajectory: 

1) The range and azimuth measurements are more accurate; and 

2) Down-linking of additional aircraft itate variables is possible 
which allows even more accurate trajectory reconstruction. 

The additional down-link capability of the Mode S transponder also 
provides the potential to improve accuracy in real-time tracking. As 
discussed in Appendix A, and in other studies [4,8], down-linking of 
aircraft state variables such as airspeed, heading, sink rate, and roll 
angle can improve the accuracy in determining position and velocity. 

The ability to determine an aircraft's position and velocity more 
precisely aids the air traffic control process in four other ways: 

1. If the controller can determine the position and velocity of the 
aircra..t he is monitoring more precisely, he can space them 
closer together. The current minimum spacing requirements 

(5 nmi enroute; 3 nmi terminal) are largely governed by the 
inaccuracies of the tracking radar and controller display. If 
aircraft can be spaced more closely, the airspace can be used 
more efficiently. (Of course, the wake vortex problem still 
must be taken into account.) 

2. By having a precise measurement of aircraft position and velocity, 
the relative dynamics between adjacent aircraft can be determined. 
This information can be used to check the possibility of close 
encounters or collision threats. This is the basis of the 
Automatic Traffic Advisory System (ATAS) which is intended to 
back up the controller. In particular, the availability of down- 
linked velocity measurements enables the establishment of 
horizontal separations, miss distance, and horizontal maneuvers 
for resolution of threat situations. Down-linked sink rate also 
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can greatly enhance the establishment of vertical separations 
and vertical maneuvers for resolving a potential conflict. Thus, 
down-linking improves airborne safety. 

3. By having down-linked velocity and sink rate, this information 
can also be used for cross-linking to adjacent aircraft. This 
information would enhance the performance of the TCAS II 
airborne collision avoidance system. The functional requirements 
and improved airborne filtering using cioss-linked information 
are the subjects of a near-future AMA study for NASA Langley. 

4. The use of down-linked information to improve ATAS or cross- 
linked information to improve TCAS enables both of these systems 
CO serve as basic sensors for Cockpit Displays of Traffic 
Information (CDTI). The advantage of down-linking to ATAS for 
CDTI purposes was studied ?nd documented in Ref. 4. The many 
applications uf the CDTI are now being determined in a joint 
NASA-FAA program. 

The above four applications plus that discussed previously for enhanced 
trajectory reconstruction for accident investigation indicate that down- 
linking of additional variables besides encoded altitude should be seriously 
explored. Fortunately, the Mode S system is designed with the flexibility 
that makes this possible. 


Ill 


MODE S TRACKING CHAR.\CTERISTICS 

Several available documents (Refs. 9-22) have been reviewer 
with respect to the Mode S Secondary Surveillance Radar (SSR) tracking 
system, its basic characteristics, its measurement errors, types of 
filtering used, anticipated utilization by the Air Traffic Control (ATC) 
system, and the data link to the aircraft. This l 'stem was formerly 
known as the Discrete Address Beacon System (DABS) . A formerly considered 
application concerning collision avoidance was known as Automatic Traffic 
Advisory and Resolution Service (ATARS). Thus, some of the recent 
literature use the DABS and ATARS acronyms. 

The ground based element of the Mode S tracking system is known as 
the Mode S sensor. Its specifications are given in Ref. 9. Further 
description is given in Refs. 10-16. The results of field testing 
three prototype versions of the sensor are found in Refs. 10,11,19-22. 
This chapter focuses on recent results presented in Refs. 19-22. 


The airborne elements of the Mode S tracking system are referred to 
as transponders. Three versions are expected to be utilized: 

a) Mode A - designed to respond to the old Air Traffic Control 
Radar Beacon System (ATCRBS) format. When the sensor signal 
is received, the transponder responds with an encoded identity 
message. 

b) Mode C - similar to Mode A, but encoded altitude is also 
included in the down-linked message. 

cj Mode S - similar to Mode C in down-linked message. However, 
when on the Mode S sensor roll call, this transponder only 
responds when its unique address is called on the sensor 
uplink message (discrete interrogation). 
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Moda A and C transponders are referred to as ATCRBS transponders. All 
three transponders will respond to both ATCRBS and Mode S sensor 
interrogations. The following sections briefly summarize what has been 
docvnnented concerning expected accuracy of the systems in terms of 
tracking a target aircraft. 

ORIGINAL PAGE IS 
OF POOR QUALITY 

Specified Airborne Transponder Errors 


The Mode S and Modes A/C (ATCRBS) transponders have built in 
specified response delays [11] as follows: 

ATCRBS: 3 + 0.5 ps (bias) + 0.06 ps (jitter) 

Mode S: 128 + 0.25 ps (bias) + 0.05 ps (jitter). 

In terms of errors, a 0.5 ps return error is equivalent to 250 ft 
range error. Thus, the above two measurement uncertainties are equivalent 
to: 

ATCRBS: + 250 ft (bias) + 30 ft (noise) 

Mode S: + 125 ft (bias) + 25 ft (noise). 

The Mode S transponder bias will remain relatively constant. The 
ATCRBS transponder bias is somewhat dependent on the received power. 

For example, biases of 188 ft, 210 ft, and 226 ft were reported as 
typical values for range of 6, 17, and 34 nmi . , respectively, [11] 
in the first flight test at the FAA Technical Center. The jitter magnitude 
(rms) remained independent of range. 

The probability of detection (Pp) witn both transponders is to be 
greater than 0.99. This is known as the blip/scan ratio (BSR) . 


Specified Ground-Based Sensor Errors 

The combination of Mode S ground sensor and Mode S airborne transponder 
is specified [9] to provide accuracy of the order 
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Range: + 150 ft (bias) + 50 ft (jitter) 

Azimuth: +0.1° (jitter) 
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Sensor only range errors measured with either Mode S or ATCRBS transponder 
reports are not to exceed + 30 ft bias and + 25 ft jitter. Thus, for a 
Mode S transponder with the specified + 125 ft range bias and + 25 ft 
noise , the above range requirements would be met . 


Flight Test Results 


Several tests were run at the FAA Technical Cevster (FAATC) to test 
the above constraints. The following sections summarize briefly the 
results of the later tests [19-22 j. 

Blip /Scan Ratio Reference 19 describes a flight test program that 
was conducted to examine the Mode S system’s ability to transfer coverage 
from sensor site to sensor site. The coverage range of the terminal area 
sensors is 60 nmi, and the enroute sensors are limited to 200 nmi range. 
Also, each sensor is limited to + 30° elevation coverage. This produces 
a "cone of silence" above each sensor. In addition, because adjacent 
sensors’ coverages do overlap, a protocol is necessary to determine 
which sensor has primary tracking responsibility for an aircraft flying 
in the range of several sites. The flight tests examined this multisite 
functioning. 

One of the outcomes of these tests was measurements of sensor blip/ 
scan ratios (BSR) . This is defined as 

„„„ Number of reports 

BSR = — - — r ^ 7 X 100. (1) 

Number of scans under track ^ ' 


Figure 1 shows a typical flight path of the test flown over tnree 
adjacent Mode S sensor sites. Terminal sites are located at the Technical 
Center, Atlantic City, and at Clementon, NJ. An enroute site is located 
at Elwood, NJ . The zenith cones are also indicated around the sensor sices. 
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Figure 1. Technical Center Surveillance Plot, Three Sensors, All 

Connected [16]. 


Table 2 is a compilation of the BSR measured at each sensor over 
several days of testing. From a total of 16,776 samples, the average 
BSR was 99.7%. This is well over the required 99.0% specification. 

Comparison Between Mode S and ARTS A second test was made at the 
Technical Center [20] to compare the surveillance performance of the 
ATCRBS mode of the Mode S sensor to that achieved with the Automated 
Radar Terminal System (ARTS) III. The ARTS III system is currently used 
in terminal areas. The surveillance performance characteristics of the 
two systems were measured using targets of opportunity. Range and accuracy 
comparisons were made by tracking an FAA test aircraft. 

The data was collected by using a common 5-foot ATCRBS antenna. The 
output of the RF portion of the sensor was alternately switched between 
the Mode S and ARTS III. 
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Table 2 


Multisite Surveillance Results [19] PAiStl IS 

CF POOK Qli!AL!TY 


Data Filter: Range 4~59 nmi; Azimuth All AZ Except 120*-145 
At Tech Center (Hangar Reflections) 


Elevation: All Elevation Angles; Test Aircraft Only (7FFFFF) 





Sample 

Size 


Blip/Scan 

Ratio 


Test 

Tech 



Tech 




No. Date 

Center 

Elwood 

Clementon 

Center 

Elwood 

Clementon 

ATC 

13 

4/2/80 

208 

261 

217 

99,5 

100. 

99.5 

100. 

24 

4/2/80 

259 

249 

177 

100. 

100. 

98.9 

100. 

15 

4/7/80 

258 

184 

123 

100. 

100. 

99.2 

100. 

18 

4/7/80 

279 

275 

1681 

100. 

100. 

95.2 

100. 

14 

4/7/80 

216 

224 

189 

100. 

100. 

100. 

100. 

26 

4/16/80 

244 


277 

99.6 


99.3 

99.6 

27 

4/16/80 

258 


262 

100. 


98.5 

100. 

19 

4/18/80 

265 

190 

193 

100. 

100. 

99.5 

100. 

20 

4/18/80 

269 

185 

197 

100. 

100. 

100. 

100. 

23 

4/18/80 

262 

N.A. 

200 

98.5 

N.A. 

100. 

100. 

30 

4/21/80 

236 



100. 



100. 

34 

4/21/80 

249 

208 

216 

98.8 

100. 

98.6 

98.8 

28 

4/22/80 

287 

240 


100. 

100. 


100. 

29 

4/22/80 

258 

213 


100. 

100. 


100. 

32 

4/22/80 


1262 



100. 


100. 

21 

4/22/80 

266 



99.6 



99.6 

25 

4/30/80 

263 

I29I 

264 

100. 

100. 

100. 

100. 

16 

5/1/80 

225 

252 

242 

100. 

100. 

100. 

100. 

21 

5/1/80 

261 

176 

285 

100. 

100. 

100. 

100. 

33 

5/1/80 


N.A. 



N.A. 


N.A. 

41 

5/5/80 

272 


301 

100. 


100. 

100. 

42 

5/5/80 

267 


296 

100. 


100. 

100. 

46 

5/15/80 

266 

116i 


100. 

100. 


100. 

17 

5/15/80 

257 

142 

232 

100. 

100. 

99.6 

100. 

52 

6/12/80 

234 


N.A. 

100. 


N.A. 

100. 

54 

6/12/80 



225 



100. 

100. 

48 

6/20/80 

288 

II2I 

283 

100. 

100. 

100. 

100. 

49 

6/20/80 

280 

N.A. 

286 

100. 

N.A. 

100. 

100, 

50 

6/20/80 

270 

271 

301 

100. 

100. 

100. 

100, 

51 

6/30/80 

280 

284 


100. 

100. 


100. 

53 

7/7/80 


280 



100. 


100. 

55 

7/7/80 

206 

240 

297 

100. 

99.6 

100. 

100. 


Total: 

7188 

4357 

5231 






Avg: 




99.8 

99.9 

99.5 

99.9 


Grand Total: 16776 

Grand Avg: 99.7 

N.A. = Read error at beginning of tape 
^ = Read error in the middle of tape 

2 = Data extraction problem 
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The Mode S sensor employs monopulse direction finding, a technique 
using a rotating fan beam antenna with a sum pattern and a difference 
pattern. The ratio of the phase and amplitudes of the signals received 
is used to determine the off-boresight angle of the target (i.e., the 
angular difference between the target position and the antenna pointing 
angle) . Reliable and improved ATCRBS surveillance data is derived from 
Mode S with a nominal 5 hits per target contrasted to today's ATCRBS 
which requires 15 or more hits per target. 

The ARTS III converts beacon video derived from the ATC Beacon 
Interrogator (ATCBI) into digital target reports. These reports go to 
the ARTS III Tracker and provide target position, velocity, altitude and 
identity to the air traffic controller. The ARTS III system uses 15 to 
18 replies per report. It also uses a sequential hit, miss counter to 
determine target detection and azimuth center marking. 

The surveillance characteristics analyzed for the two systems 
include BSR, identity (Mode 3/A code) reliability, altitude reliability, 
and extra (false or split) targets. The systems were evaluated with both 
dissemination and no dissemination of uncorrelated beacon data. Only 
correlated data results are reported here. Targets were limited from 4 
to 45 nmi in range, 0.5® to 30® in elevation, and those azimuths where 
no ground reflection interference exists. 

The results of 24 surveillance tests are given in Table 2. As 
seen from the mean values for the 24 data sets, the Mode S (ATCRBS Mode) 
is always superior to that of ARTS III. In particular. Mode S altitude 
reliability is 6% better than ARTS III. 

Range and azimuth accuracy of the Mode S (ATCRBS mode) and ARTS III 
were determined using Technical Center aircraft flying predetermined 
radial and orbital flightpaths while being tracked by the Technical Center 
Nike-Hercules precision tracking facility. The range and azimuth positions 
of the test aircraft, as detected by both the Mode S and the ARTS III, 
were compared to the Nike-Hercules positions to establish accuracy 
measures. Both inbound and outbound radial flights and clockwise and 
counterclockwise orbital flights were used to collect the data. 
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Table 3. Mode S ( ATCRBS Mode) Correlated Beacon Reports Compared with ARTS III Reports [17] 
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Value 98.7 96.3 99.6 97.8 96.2 90.3 0.68 1.65 

26 Data 


OF POOR QUALnY 

Table 4 presents the resulting azimuth and range test errors of 
Mode S compared to ARTS III. The Mode S mean a«:imuth error is an order 
of magnitude better than ARTS III, and its standard deviation is 5 to 
20 times better than ARTS III. This is because (a) Mode S uses a 14 bit 
antenna shaft encoder (precision of 0.022®) rather than the ARTS III 
12 bit encoder (precision of 0.088®), and (b) Mode S uses the monopulse, 
of f-boresight angle correction technique rather than the ARTS III 
sequential hit-miss counter. 


The least significant bit of the ARTS III range is 380 ft, so the 
expected deviation in range due to quantization effects is 110 ft. The 
Nike-Hercules tracker used as a reference has a 10 ft random error. 

From Table 4, the ARTS III range error varied from 128 to 136 ft. 


Table 4. Comparison of Mode S (ATCRBS Mode) and ARTS III Accuracy [20] 


•Azimuth Errors (Degrees) 



Mode 

S (ATCRBS 

Mode) 


ARTS III 



Sample 


St andard 

Sample 


St andard 


Size 

Mean 

Deviation 

Size 

Mean 

Deviation 

Outbound Radial 

78 

0.003 

0.041 

78 

0.049 

0.198 

Inbound Radial 

137 

-0.004 

0.029 

121 

0.072 

0.255 

Clockwise Orbit 
Counterclockwise 

147 

-0.004 

0.025 

149 

0.117 

0.238 

Orbit 

41 

0.028 

0.018 

39 

-0.098 

0.312 




Range 

Errors (Feet) 





Mode 

S (ATCRBS 

Mode) 


ARTS III 



Sampl e 


Standard 

Sample 


Standard 


Size 

Mean 

Deviation 

Size 

Mean 

Deviation 

Outbound Radial 

78 

69 

27 

78 

11 

128 

Inbound Radial 

137 

87 

33 

121 

-126 

121 

Clockwise Orbit 
Counterclockwise 

147 

67 

27 

149 

-4 

134 

Orbit 

41 

70 

26 

39 

-70 

136 
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The least significant bit for the Mode S range is 60 ft causing the 
expected deviation due to quantization to be 17.3 ft. From Table 
the Mode S range error varied from 26 to 33 ft. 

The data from Table 4 is presented in graphic form in Figs. 2 
and 3. In summary, the mean and standard deviations in comparing 
these systems were as follows : P'Vl'J ^3 

OF POOR QUALITY 

Range (ft) Azimuth (deg) 




Mean 

Std. 

Mean 

Std. 

Mode 

S 

+ 74 

29 

+ .006 

0.029 

ARTS 

III 

-46 

129 

+ .035 

0.243 


The significant imp-overaent of the Mode S sensor in tracking ATCRBS 
trarsponder equipped aircraft is apparent. 

Comparison of Mode S Sensor Accuracy with Mode S and ATCRBS Transponders 
The purpose of the third set of tests was to evaluate the Mode S sensor 
accuracies at the Technical Center, Cleme ..;on, and Elwood, NJ , sites [21]. 

In all cases, the test aircraft was instrumented with both ATCRBS and 
Mode S transponders. 

The monopulse azimuth accuracy was verified by using special test 
transponders at each sensor site. These transponders, referred to as 
calibration performance monitoring equipment (CPME), are permanently 
installed at surv'^syed locations within the coverage area of the associated 
sensor. Azimuth calibration bias adjustment is accomplished by reference 
to surveyed azimuth value for the CPME location. 

Figure 4 illustrates the horizontal profiles of the aircraft 
flights. The radial flightpaths were flown at azimuth angles that align 
the Nike tracker and the Mode S sensor. For the orbital flights, nearly 
constant slant range was maintained around the sensors. Orbital altitudes 
were varied so that elevation angle ranged from 2° to 30°. 
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Figure 2. Mode S (ATCRBS Mode) ARTS III Azimuth Accuracy Comparison [20]. 





MODE S/C«TS ACCURACY COMPARISONS AIRCRAPT N-EO 



Figure 3. Mode S (ATCRBS Mode) ARTS III Range Accuracy Comparison [20]. 








The specific results for the Technical Center sensor are presented 
here. These are similar to the data results of the other terminal sites. 
Table 5 summarizes the azimuth and slant range residuals for Mode S 
and ATCRBS transponders on each of three flight tests. For each date, 
the total number of samples, mean residual, and standard deviation were 
tabulated. 

Figures 5 and 6 show the two transponder test result histograms 
for the first flight date. These histograms illustrate typical raw 
data distribution and extreme values. Most of the histograms tend to be 
bell-shaped except for Mode S azimuth residuals. The Mode S azimuth 
residuals are skewed to the left because of a reduction in accuracy at 
the higher elevation angles due to the antenna beam pattern and the Mode S 
technique of requiring only s. single reply per scan. The ATCRBS azimuth 
residual histograms show only minor skewness because it averages the two 
replies closest to the boresight. 

Typical azimuth residuals as a function of elevation angle are 
shown in Fig. 7 for both Mode S and ATCRBS transponders. Figure 8 
shows range residuals as a function of slant range. A range-dependent 
bias is also seen here. 

Emperical mathematical models were developed for azimuth residuals 
based on the secant of the elevation angle and range residuals based on 
slant range. Table 6 tabulates the results of the azimuth and range 
residual regression models for the Technical Center sensor. These model 
results for the first day are also illustrated in Figs. 7 and 8. 

The Mode S azimuth model in Fig. 7 reduces the standard deviation from 
0.062° to 0.038°. For range, the ATCRBS data shows about a 1.5 ft change 
in range residual per nmi. 

One problem with the sensors is that they exhibited both long and short 
term drift in the azimuth residual data. Figure 9 illustrates azimuth 
residual change for the Technical Center sensor. The CPME data is super- 
imposed. The CPME data display the same general shifting of bias as the 


Table 5, Technical Ce iter Sensor Azimuth and Slant Range Residuals [21] 
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7/1 19,000-ft radials 10 1,299 -39 30 1,289 108 
7/23 39,000-ft radials 8 635 -79 42 624 83 
7/22 Orbital arcs 13 1,043 -114 22 1,033 33 
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Figure 6. Azimuth Residual Data Histograms for Technical 

Center Sensor Radial Flight of July 1, 1980 [21]. 
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Table 6. Mathematical Regression Models of Azimuth and Range Residuals [21] 


ORtGIWAL PAQS 
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Figure 8. 


Range Residual Mathematical Model Plots 
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Figure 9. 


Mode S and ATCRBS Azimuth Residuals Illustrating Long 
and Short Term Azimuth Drifting [21]. 
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azimuth residual data; this indicates that the problem is inherent in the 
sensor. Figure 10 illustrates the azimuth residual results for the two 
terminal sensors for elevation angles less than 12.5°. 

Comparison of Mode S and ARTS for Range and Azimuth Resolution The 
latest flight tests [22] were conducted at the Technical Center to compare 
range and azimuth resolution capabilities of a Mode S sensor with ARTS III 
for the following scenario: 

a) Two aircraft fly close together in both range and azimuth; and 

b) Both aircraft are equipped with ATCRBS transponders. 

The minimum achievable range and azimuth separations of two aircraft 
without garbling of either aircraft's identity (A) code were determined 
for both systems. The resolution results were compared to positional 
aircraft separation data, collected concurrently by the precision Nike- 
Hercules tracking radar to determine the relationship between A-code 
garbling and aircraft separation. 

Note here that two targets with Mode S transponders can always be 
separated because of their .lique addresses. Similarly, two targets - 
one with Mode S and the other with ATCRBS transponder - can be resolved 
because they are interrogated at different time periods. So, garbling 
is only a problem with two ATCRBS equipped targets. Garbling is caused 
when two RF signals, while propagating through space, overlap or interfere 
with each other. This occurs when two aircraft are so close that there 
transponders reply to the same ground interrogation at nearly the same time. 

Figure 11 illustrates the radial flightpaths of the two test air- 
craft. The target flew the 110° radial, and the overtaking aircraft flew 
at either 110°, 111°, or 112°. They were separated in altitude by 500 ft, 
and the separation ranged from 0° to 4° in azimuth and 0 to 20,000 ft in 
slant range. 

Figures 11 and 13 illustrate the resolution results for the Mode S 
(correlated-only reports) and ARTS III sensors. Each cell in the matrix 
contains the total number of samples available and the percentage of 
resolved reports. 
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ZIMUTH RESIOURLS FOR 



Technical Center and Clementon Sensor; 














The flight test results iudicato that an 89°: beacon resolution was 
achieved with the Kode s sensor in the aircraft separation intervals of 0° 
to 2° in azimuth and 0 to 10,000 ft in ranpQ. For ARTS III, the same 
aircraft separation intervals were b2;'o resolved. The minimum achievable 
range separation, without garbling, was about 10,000 ft for both systems. 
The azimuth separation was 2° for the Mode S system and 3.2° for the 
ARTS III system. 



Figure 11. Radial Flightpaths of each Aircraft [22], 
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ENCODING ALTIMETER CHARACTERISTICS 


The barometric altimeter determines altitude based on measured 
static pressure. If the ground pressure setting on the altimeter is set 
at 29.92 in Hg (standard sea level pressure), the altitude reading is 
referred to as pressure altitude. This altitude is encoded to the nearest 
100 ft level and linked to the ground with the Mode C transponder. 

There are several sources of error for this measured altitude. 

These can be grouped into two categories : 

1. Meteorological error due to the atmosphere being non-standard; and 

2. Altimeter instrumentation errors. 

Both of these errors are discussed in turn. 


Meteorological Error 


The baro-altimeter is calibrated to convert static pressure to altitude 
based on the Standard Atmosphere. The Standard Atmosphere [23] is a model 
of pressure, density, and temperature as functions of geopotential altitude. 
This model is based on the following assumptions: 


1. Sea-level conditions are as follows: 

Pressure p = 29.921 in Hg = 2116 Ib/ft^ 

o 

= 101325 N/m^ = 1013.25 mbars . 

3 

Density = 0.0023769 slug/ft 


Temperature T^ = 59.0°F = 288.16^K = 15.0°C 

2 

Gravity g^ = 32.17 ft/sec (assumed constant for 

geopotential altitude) 




ir»v 4 7. 
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2 . 


The equation of state for dry air holds: 


PASS IJ5 
OF POOK QUALITY 


p/p = 8 q RT (2) 

where R is the gas constant (R=53.3 for air) and T is in 

‘’Rankine (°F + 459.4'’), Also, p and p are pressure and 
2 3 

density in Ib/ft and slug/ft , respectively. 

3. Temperature varies linearly with altitude up to the tropopause 
(at 36,089 ft = 11 km) . 

T (“F) = 59 - .003566“h (3) 

Thereafter, it is constant at -69.7‘’F up to 65,617 ft (20 km). 
Equation (3) is written as 

T = T^ - ah , (4) 

where a is known as the lapse rate. 

4. The differential pressure over a small slice of air can be 

written as 

dp dh 

p " ” RT ’ 

or 

dp = -p dh , 

where dh is the width of the slice. 


(5a) 

(5b) 


Using Eqs. (4) and (5) gives 


ff ■ / 


dT 
a RT 


or 


In 


(t) ■ (^) ’ 


From Eqs. (4) and (7), we can write 


(P/P„)^“ - (I„ - ah)/T^ , 


( 6 ) 

(7) 

( 8 ) 
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and 



Equation (9) is the expression of geopotential altitude as a 
function of lapse rate a, sea-level conditions p^ and and measured 
pressure p. It is the basic equation of the baro-altimeter. 

Altimeters are designed so that they can be corrected for deviations 
in surface pressure. The pressure is measured on the ground, and an 
equivalent sea-level pressure p^^ is computed so that Eq. (9) holds for 
the known altitude h of the measurement site (usually the runway) and the 
measured pressure p. This corrected pressure p^^ is entered into in the 
altimeter setting whenever the aircraft flies below 18000 ft. This term 
is then used by the altimeter so that it reads runway altitude upon 
landing. Above 18000 ft, the standard p^ (29.92 in Hg) is used. There 
are no means for correcting for deviations in the lapse rate a or ground 
temperature T^ in most altimeters. Thus, these variations cause an error 
between true altitude and the indicated baro-altitude. 

Figure 14 compares standard atmospheric pressure and temperature 
vs geopotentiai altitude with the average values measured each day in 
Albuquerque as an example during September 1980. The deviation in 
temperature is most pronounced. However, if a new lapse rate a, reference 
temperature T^, and pressure p^ are defined for each segment, then Eqs . 

(4) and (7) hold for this period. The vertical distance between 
measured and standard pressure lines result in the average altitude error 
as a function of altitude. The maximum vaxne of this error is about 600 m 
(2000 ft) at 10 km altitude. 

The arrows in Fig. 15 show the lowest and highest mean monthly tempera- 
tures obtained for any location between the equator and pole. Estimates of 
the one-percent maximum and minimum temperatures that occur during the 
warmest and coldest months, respectively, in the most extreme locations are 
shovm by dashed lines. Values below 30 km are based on radiosonde observations. 
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Fi^’,ure 15. Range of Systematic Variability of Temperature Around 
the U.S. Standard Atmosphere, 1976 [23j. 


Values shown for the various levels by envelope curves could not 
possibly be encountered at all altitudes at a given location and time. The 
warmest layers near the surface, for example, are associated with the 
coldest temperatures at the tropopause. 


At latitudes between 30 and 90®N, maximum mean monthlv temperatures 
at altitudes below 25 Itm usually occur in June or July, and minimum values 
occur in December or January. 


If one can obtain the temperature profile and the surface pressure, 

then Eq . (9) can be used to compute the actual geopotential altitude for 

a given day from measured pressure p. This can be used with Eq . (';) , 

based on standard dav values of (a, p , T ), to determine the altitude 

o o 

error due to meteorological conditions. 

For higher altitudes, the assumption inherent in Eq . (5b) that the 

gravitational acceleration remains constant at the sea-level value begins 
to produce another minor source of error. Note that the altitude h of 
Eqs. (5) is referred to as "geopotential altitude". The Newtonian 
equation for gravitational potential is 
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( 10 ) 



where r is the earth's mean radius 
e 

the true, or "geometric altitude", 
atmospheric balance is 


(r = 3959 s.mi = 6371 km), and z is 
The actual hydrostatic equation for 


dp = -pg dz, 

which replaces Eq . (5b) . From Eqs, (5b), (10), and (11), we get 

the relationship that 


( 11 ) 


( 12 ) 


At the tropopause (h = 11 km), the difference in h and z is 19 m (62 ft). 
The geopotential altitude is used with the definition of the Standard 
Atmosphere as a matter of convenience. 


Currently, the FAA Technical Center is running tests where the true 
altitude is being measured and compared to the encoded altitude linked to 
the ARTS III tracking facilities [24]. This is being done by tracking 
enroute aircraft with the Wallops Island precision tracking radar and comparing 
the derived altitudes with those reported at Washington, Baltimore, and 
Philadelphia ARTS facilities. The results are to be published in May 1983. 
These results will produce more examples of iaherant meteorological 
altitude errors. 


Altimeter Instrumentation Errors 


Altimeter instrumentation errors can be classified into three 
categories : 


1. Static pressure errors (.also referred to as static defect errors); 

a. Dynamic pressure addition caused by speed and angle-of- 
attack being applied to the pitot tube; 

b. Trapped condensation or foreign material blocking the 
pneumatic system; 
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c. Leaks in the static pressure system; 

d. Deformation to pressure source or variations in 
mounting angle and aircraft skin. 

2. Mechanical errors of the instrument; 

a. Diaphrani or scale error due to the nonlinearity in the 
aneroid and linkage system; 

b. Instability - tendency of original calibration to change 
after repeated ascents and descents; 

c. Friction in mechanism; 

d. Hysteresis due to imperfect elastic diaphragm; 

e. Gear backlash; 

f . Temperature variation; 

g. Lag in transmission of pressure in tubing when in climb 
or descent. 

3. Reporting eriors; 

a. Encoding aparatus error; 

b . Quantization error. 

These three error sources are depicted in Fig. 16 with slightly 
different sub-categories [28] . It is very difficult to compute how 
each of these types of errors contribute to the overall error of an 
individual altimeter. However, there have been many tests run [25-28] 
to determine the nature of altimeter errors, and other tests are in process 
[27]. The following is a summary of material on altimeter instrumentation 
error characteristics as documented in Refs, 25-28. 

Reference 26 compares the altitudes read for a pacer aircraft 
(calibrated by a phototheodolite) with that of 115 subject aircraft, each 
flying (one at a time) beside the pacer aircraft. These were general 
aviation aircraft flying at 4500 and 8500 ft. All aircraft were flown 
with altimeters set at 29.92 in Hg. Thus, all altimeters were reading 
pressure altitude. Both error in indicated altitude (displayed by the 
barometric altimeter in the cockpit) and reported altitude (used by ATC on 
the ground) were measured. 
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Figure 16. Components of Aircraft Altimetry Instrumentation [28]. 






of the 115 aircraft tested, the overall error in indicated altitude had 
a mean of -3 ft and a standard deviation of b8 ft. This is depicted as a 
histogram in Fig. 17. From the aircraft tested, A5 were calibrated 
regarding their altitude reporting error. This error had a mean value of 
33 ft and a standard deviation of 111 ft; it is shown in Fig. 18. 


For air carriers, the following range of standard deviations in 
altitude error were reported [2b]. 


Indicated errors: 

ATCAC Report : 

ICAO Report: 


o, = 31 ft 
A 

a, = 43 - 54 ft 
A 
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Installation errors: 

FAR requirements: 62.5 ft 

SRI Report: = 56 ft. 

Here, a is the installation errors in indicated altitude (referred to as 
c 

correspondence error"). The root sura square with a of 54 ft and a of 
5b ft gives a total of 78 ft. 


None of this information gives an indication of the dynamic characteristics 
of the errors. However, Billroann [27] devised a model for general aviation 
altimeter errors based on data taken from two flights at 2500 ft. The 
sample mean, variance, and standard deviation of encoded altitude error 
were 0.3 ft, 579.6 ft" and 24.1 ft, respectively. These errors proved to 
have close to a normal distribution. The measurements were sampled every 
one sec and were highly correlated. The expression describing the correlated 
altitude errors E was 

E(t) = E(t-l) + E(t-2) + W; t ^ 2, (13) 

2 

where W is zero-mean white noise with a variance of 111.1 ft", and and 4>2 
were 1.066 and -0.191. Again, the standard deviation of E(t) was 24.1 ft. 

If the altitude is sampled every 5 sec, this correlation can be ignored, 
and the error can be treated as white noise. 
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from flight test data 



ERROR IN FT. 


Figure 17. Average Errors in Indicated Altitude for the GA 

Population [26 J. 




FREQUENCY 




pg peOll QUALITY 


NOTE: DATA OBTAINED FROM FLIGHT TESTS 

SAMPLE SIZE » 45 FT. 

MEAN (^i) = 33 FT. 

STANDARD DEVIATION (cr) - 111 FT. 


ERROR IN FT. 


Figure 18. Histogram of Average Error in Reported Altitude for 
the General Aviation Population [26]. 



In Ref. 22, several researcners explain experiences with extensive 
flight tests of altimeters in military aircraft. In particular, the 
altimetry errors due to variations in the pitot static system were 
assessed because this was judged to be the source of a large uncertainf; . 
The factors of the pitot static tube that affect the indicated altitude 
are those listed as "static pressure errors" noted previously. Ninety- 
four different military aircraft were tested for a total of over 1650 
flight hours. Aircraft tested included the F-104, F-4C/D, T-38, and A-7D, 
of which there are 1800 aircraft. 

Figure 19 shows the variation in altitude corrections as a function 
of Mach number for four F-104' s. As much as 200 ft difference is seen 
between aircraft altitude readings at high Mach. 

Figure 20 shows altitude correction variations for four F-4C's. 
Differences in altitude of more than 300 ft are seen. In fact, 10 % of 
the static ports inspected on 25 F-4C/D aircraft did not meet technical 
order tolerances; this produced the indicated bias variations in the 
altimeter readings. 

Figure 21 shows the effect of pitot tube burrs on a T-38 aircraft 
altitude reading (180 ft error) . The author concluded that the chances of 
a T-38 aircraft meeting a + 250 ft accuracy requirement with its altimeter 
were small. 

Figure 22 shows A-7 aircraft altitude reading variations; as much as 
300 ft altitude variation is seen. 

From these examples, it is seen that large speed dependent biases can 
exist on the altimeter output. Calibration checks of the pitot tubing are 
continuously necessary to minimize this error source. 

The ARINC 575 performance Requirements for altimeter accuracy are 
shown in Fig. 23, along with results of flight tests of a commercial 747 
125J. The point is that altimeter-indicated accuracy can be kept small 
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NOTE: Position error correction curves at 36,000 feet. 


Figure 19. F-104 Pilot Static System Data [25]. 
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Indicated Mach Number 


Figure 20. F-A Flight Test Data [25]. 
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Test Data [25] 




Deviation* (F 



Fij’ure 23 . Flight Test Altitude Deviations for Aircraft Carriers [ 25 ] 
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(x^^ichin the ARINC requirements), but calibration maintenance of the pitot 
tubes is required. 

Another study of a 7A7 altimeter was made [25] to determine typical 
dynamic lag characteristics due to descending or climbing flight. The 
lag for the altimeters grew from about 1 sec at 10,000 ft altitude to 
about 3 sec at 35,000 ft. This deviation can be compensated. For a 
simulated emergency descent at 17,000 ft at 170 ft/sec, the altitude error 
was 270 ft. After lag correction, this was reduced to 30 ft. 

As a final note, we must consider the effect of quantizing the 
reported altitude to the nearest 100 ft accuracy. The reported altitude 
goes through the following sequence of manipulations: 

a. The measured altitude is subject to certain mechanical 
errors, as previously discussed. 

b. encoded, based on the standard 29.92 in Hg setting. 

This report is quantized f, o the nearest 100 ft. 

c. The ARTCC computer applies a local pressure correction to 

to convert it to pressure altitude. The result is 
quantized again to z^. 

The reported altitude error is the difference between z and z^ , the true 

c t 

pressure altitude. The quantization error in both z^ and z^ has a standard 
deviation of 50//J”. These are each root-sum-squared with the pre-quantizer 
error (which is similar to the indicated altitude error). For the data 
taken in Ref. 26, (see Fig. 18) the pre-quantizer error had a standard 
deviation of about 100 ft. Reference 28 explains that there are both 
encoding and quantization errors in the encoding process. These are 

depicted in Fig. 24. For •<...■ carriers, the encoding error has a limit 
of 15 ft (2a). For geneia* aviation, FAR Part 91.36 requires that the 
quantized altitude, as used in the Mode C altitude report, correspond to 
the indicated altitude to within + 125 ft (2a) when installed. 
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Reference 28 categorizes the instrumentation errors shown in Fig. 

16 according to whether they are tor air carriers or general aviation. 

The air carrier altimetry equipment is l.argely conti oiled by ARINC 
characteristics (specifically, ARINC characteristic 545) , while the 
general aviation fleet is largely controlled by t'le Federal Aviation 
Regulations (specifically, Part 43 Appendix E) . For the air carriers, 
the static source errors were examined in terms of Mach, angle-of-attack , 
calibrati'.n , and aircraft variability effects. The air carrier error 
estimates, in terms of standard deviations, are shown in Table 7. 

Static error is the largest contributer to general aviation altimeter 
instrumentation error. These errors are reduced if the altimeter has 
static defect correction (SDC) , although this feature is relatively 
expensive. The error estimates for general aviation are shown in Table 8. 

In conclusion, there are various references concerning modeling of 
altimeter instrumentation error magnitudes. Many factors influence the 
reported error magnitudes including manufacturing specifications, care in 
installation, type of aircraft, flight conditions, and regularity of 
calibration. For this reason, the analyst has a large range of choices in 
deciding what magnitudes to assign to various altj.neter instrumentation 
errors . 

For simulation purposes, the altimeter can be modeled as having 
additive white noise of zero mean and standard deviation ranging between 
25 ft and 200 ft. To this is added a scale factor error proportional 
to altitude and speed. At 30,000 ft this could be up to 300 ft at Mach 
0.8. In addition, there is an altitude dependent lag term (up to 3 sec 
time constant) affecting climb and descent readings. Finally, to these 
errors, one must add two quantization error sources. The choices depend 
on whether the user is concerned with general aviation or air carriers, 
and whether a conservative or optimistic accuracy model is desired. 
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Figure 24. Mode 0 Encoding and Quantization Errors [28] 


Table 7 
System 

. Estimated 
Performance 

Standard Deviation 
Among Air Carrier 
U.S. Fleet [28], 

in Total 
Ai rcraf t 

Altimetry 
of the 

ALTITUDE(MSL) STATIC 

INSTRUMENT 

MODE C 

TOTAL RSS (feet) 

(feet) 

SOURCE 



w/o Mode C 

v/ Mode C 

BL 

31 

12 

30 

33 

45 

5 I 

36 

12 

30 

38 

48 

10 t 

41 

12 

30 

43 

52 

13 t 

46 

19 

30 

50 

58 

20 t 

52 

25 

30 

58 

65 

25 t 

57 

31 

30 

65 

71 

30 t 

61 

38 

30 

72 

78 

35 t 

66 

44 

30 

61 

66 

40 t 

75 

50 

30 

90 

95 
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Table 8. Estimated Standard Deviation in Total Altimetry 
System Performance Among General Aviation Aircraft of the 

U.S. Fleet [28], 


Ai,TITUDE(MSL) 

(feet) 

STATIC 

tRSioP 

XRSTRUHEKT 

QUAi^Tl- 

NATION 

TOTAL RSS 
v/o Mode C 

is feet 
w/ Mode C 

SL 

78 

26 

63 

82 

104 

5 1 

95 

29 

63 

99 

118 

XO K 

X09 

38 

63 

115 

132 

15 K 

125 

45 

63 

132 

147 

20 t 

140 

52 

63 

149 

162 

25 t 

X55 

61 


166 

178 

30 t 

168 

69 

63 

182 

192 

35 t 

185 

77 

63 

200 

210 

40 K 

200 

86 

63 

216 

227 
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ERROR MODEL DE ELOPMENT 


Chapters III and IV documented literature findings concerning the 
accuracies of measurement of (a) slant range and azimuth from the new 
Mode S tracking system, and (b) altitude as measured by a baro-altimeter 
and encoded for transmission via the airborne transponder. The purpose 
of this Chapter is to summarize these findings in terms of measurement 
error models. The rationale is given for choice of specific error 
characteristics. 


Slant Range 

The ground tracking system, referred to as the **sensor”, may be 
either the Mode S design or the current ATCRBS design (ARTS III in the 
terminal area; ARTCC enroute) . The airborne transponder also may 
be of the Mode S design or the current ATCRBS design (either Mode A or 
Mode C) . From these combinations, three different possibilities, or 
configurations, exist with regard to tracking accuracy: 

1. Mode S sensor and Mode S transponder; 

2. Mode S sensor and ATCRBS transponder; and 

3. ATCRBS sensor and Mode S or ATCRBS transponder. 

We refer to these in the following as Set 1, 2, or 3, respectively. 

It is first necessary to examine the possibility that a target return 
is received if it is scanned. This blip/scan ratio (BSR) for the three sets 
was experimentally measured to be 
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Set 

BSR 

HI 


1 

99.7 

(Table 

2), 

2 

98.7 

(Table 

3), 

3 

96.3 

(Table 

3). 



There is some uncertainty in these values because of limits to the amount 
of data. However, simulation we can approximate these Set values 

as 99.5%, 99.0%, and 96.5%, respective!’;^ (or any other set we wish to 
choose). 

The BSR factor can be mathematically modeled as uniformly distributed 
random numbers between 1 and 1000. If the random number is greater than 
965, say, for Set 3, then no report is received. 

The second factor is the reliability of the return. That is, a return 
may be received, but dropped bits or noise may cause it to not be correlated 
with similar measurements in a tracking sequence. This would cause it to 
be rejected. We assume that this is similar to Mode A (Identity Code) 
Reliability as reported in Table 3. Thus, 

Set Reliability 

1,2 99.6 

3 97.8 

Again, this reliability factor is modeled as a unifoirm distribution. The 
BSR and reliability factors govern range and azimuth measurements, alike. 

Another factor to be noted is the proximity of the tracked aircraft 
with respect to any adjacent aircraft. If two proximate aircraft ha^e 
ATCRBS transponders, there can be tracking losses if they are too close 
(2° - 3.2° in azimuth and 10,000 ft in range). The Mode S sensor is 
superior to ARTS III for this particular problem, as described in Chapter III. 
Usually, this does not have to be considered, because we are primarily 
interested here in the tracking accuracy of a single aircraft. 
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With regard to range accuracy, the next factor to consider is the 
position of the tracked aircraft with respect to the sensor site. Terminal 
sensors have 60 nmi range, and enroute sensors have 200 nmi range. The 
elevation upper limit is 30°, and the lower limit is 0.5° above the 
horizon. Of course, buildings and other terrain features must be taken 
into account in defining the horizon. 


The range error has three components : 

1. Independent white noise; 

2. Range dependent bias; and 

3. Quantization effect. 


From the test data discribed in Chapter III, the following overall 
accuracies were reported: 

Set Mean (ft) Standard Deviation (ft) 


-39 

30 

(Table 

5) 

108 

44 

(Table 

5) 

74 

29 

(Table 

4) 

-46 

129 

(Table 

4). 


The Table 5 values here are only one case. For the first two cases 
above, mean range-dependent (bias) errors were fit to the data with the 
resulting bias being: 


Set 

1: 

b 

r 

—66 + 0 . 83 * ^SL * 

(14) 

Set 

2: 

b 

r 

61 + 1.48 * , 

(15) 


where slant range r„. is measured in nmi. This correction lowered the slant 

O Li 

range standard deviations to 27 ft and 38 ft, respectively. 


If this range dependent bias is known at a particular site, it is 
worthwhile to account for this bias. Otherwise, the raw range measurement 
errors can be added as a constant bias b^ plus white noise n^: 


"SLin 


= r 


SL 


+ b + n 


(16) 
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Here, b is a constant bias with value + 100 ft, and n is white noise 
r — ’ r 

with typical standard deviations of 40 ft for Sets 1 and 2 and 130 ft for 
Set 3. 


The simulation of individual range measurements can also account for 
the quantization effects. For Set 3, the range measurement is quantized to 
the nearest 380 ft. This produces an equivalent standard deviation of 
110 ft. Thus, for Set 3 with a total standard deviation of 130 ft, say, 
the quantization would contribute 110 ft, and the white noise standard 
deviation would be 69 ft. 

For Sets 1 and 2 (Mode S sensor) , the range quantization is 60 ft with 
an equivalent standard deviation of 17 ft. Thus, for 30 ft total standard 
deviation, 17 ft would be from quantization, and 25 ft would be the white 
noise standard deviation. 


Azimuth 


Azimuth measurements also can be classified in the same three sets as 
range. Also, the BSR, return reliability, proximity of another aircraft, 
and coverage volume factors are the same as for range. In addition, the 
azimuth error has the three components: 

1. independent white noise; 

2. Elevation dependent bias; and 

3. Quantization effect. 

From Chapter III, the following accuracies were listed: 


Set 

Mean (deg) 

Standard 

Deviation 

1 

-0.032 

0.062 

(Table 5) 

2 

-KJ.036 

0.030 

(Table 5) 

2 

+ .006 

0.029 

(Table 4) 

3 

+ .035 

0.243 

(Table 4) 
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For the first two cases above, mean elevation dependent (bias) errors were 
fit to the data with the resulting biasjjeing: 


Set 1: 
Set 2: 


= 1.58 - 1.59 secant (Eil) , 

= 0.27 - 0.23 secant (Ei). 


(18) 


These curve fits lowered the azimuth standard deviations from 0.062® to 
0.038® for Set 1 and 0,036® to 0.030® for Set 2, Again, if this effect 
is known for a particular site, it is wise to account for it. Otherwise, 
the azimuth measurements can be modeled as 

a = a + b + n . (19) 

zm z a a 

Here, b is a constant bias with value of + 0,04°, and n is white noise 
a ■*“ a 

with standard deviations of 0.05° for Sets 1 and 2 and 0.25° for Set 3. 


The Mode S sensor azimuth measurements are quantized to 0.022° which 
contributes 0.006° to the noise error. The ARTS III sensor is quantized 
to 0.088° which contributes 0.025° to the noise error. 


Encoded Altitude 


Here we are modeling the difference between the aircraft's true 
altitude above sea level and that reported by the encoding altimeter through 
the transponder. 

First, consider the meteorological errors. If T^^, and p^^ are 

reference temperature, lapse rate, and pressure on a given day, the 
geopotential altitude will be 



This is based on Eq. 


(9). Thus, a bias type error will be of the form 



( 21 ) 
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The difference in geometric (true) altitude and geopotential altitude 
h is found from Eq. (12) to be • 


Ah2 



( 22 ) 


The measured altitude also has white noise n, , scale factor error 

h 

proportional to altitude and speed (e^ + e^M) and the quantization effect. 
The altimeter reading would be 


h^ = (1 + + £2^) h - Ah^ - Ah2 + n^ (23) 

This would be quantized to the nearest 100 ft. A correction term to 
account for pressure error at the sensor site would be added, and the 
result would again be quantized to h^^. In Eq. (23), n^ would have a 
standard deviation of 50 ft for example, and e^, would be set to produce 
up to 300 ft error at Mach 0.8 and 30000 ft. 


If the aircraft is climbing or descending, Eq. (23) would be 
modified to account for pilot tube lag effects. 


h 


mo 



h ) . 
mo 


(24) 


Here, h^^ is the output to the quantizer and is the time constant of 
the tubes (1-3 sec proportional to altitude). From Eq . (24), we can write 


h 

mo 


h 

m 

T, s+1 • 

n 


(25) 


where s is the Laplace operator. 


The signal h is received on the ground. To this, a local altitude 
mo 

correction term is added to account for variations from standard day 
atmosphere. The result is again quantized to the nearest 100 ft. 
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SUMMARY AND RECOMMENDATIONS 
Summary 

In Chapter II, material was presented from the FAA 20-year plan 
which indicates that the U.S. airfleet is projected to continue to grow 
to over 400,000 aircraft by 2000. Currently, about 60% of the GA aircraft 
are transponder equipped which allows good tracking by the ATCRBS systea. 
This percentage is expected to grow, and there are strong motivating 
factors to encourage transition to Mode S transponders. These include: 

1. It may be mandatory to have a Mode S transponder to 
receive ATC clearance above 6000 ft by year 2000. 

2. The Mode S transponder will enable the user to receive many 
new services via the data link which are now under development. 

Thus, there are both the "carrot and stick" factors pushing for acceptance 
and growth in use of Mode S. 

The final section of Chapter II summarizes four additional reasons 
why down-linking of aircraft velocity and sink rate (in addition to the 
current encoded altitude) is a promising idea. These reasons, plus the 
application for enhanced trajectory reconstruction for accident investiga- 
tion, indicate that down- linking of additional variables should be 
seriously explored. 

From Chapter II, the Mode S system is seen to perform better in 
flight test than required by its specifications. The blip-scan ratio 
was 99.7%. Signal reliability was 99.6%. These compare to 96.3% and 
97.8%, respectively, for the ATCRBS system. The Mode S azimuth measure- 
ments have accuracies with noise standard deviations of less than 
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0.05'" compared to 0.25® for ATCRBS. The Mode S azimuth measurement has 
a bias error which is elevation dependent and should be taken into account. 
The Mode S range measurement has a standard deviation of about 40 ft 
compared to 130 ft for ATCRBS. 

From Chapter IV, the encoded altimeter measurements have two types 
of errors: (1) Meteorological error due to the atmosphere being non- 
standard, and (2) Instrumentation errors. The meterological error can 
be greater than 2000 ft at higher altitudes, but it can be compensated 
by knowing the temperature vs altitude profile and surface temperature 
and pressure. The instrumentation errors include (a) static pressure 
errors, (b) instrument mechanical errors, and (c) encoding errors. These 
errors can be modeled as noise with standard deviation in the range 25-200 
ft, and an altitude and speed-dependent scale factor error causing up to 
300 ft error. The altimeter errors have a wide range of magnitudes, and 
they depend on type of aircraft as well as altitude and speed. 

Chapter V presents mathematical models of the ATC tracker and 
encoding altimeter errors. These are suitable for simulation purposes. 

In Appendix A, analytical methods were used to compare Mode S tracking 
accuracies in the horizontal plane to those currently available with ATCRBS 
tracking. Two error sources were investigated - (a) the azimuth measurement 
error modeled as white noise, and (b) the deterministic acceleration error 
caused by the tracked aircraft turning. Also, the effects of range to the 
target were included, and range was varied from 15 nmi to 60 nmi. The 
standard deviation of the azimuth error was 0.04° for Mode S and 0.25° for 
ATCRBS. The tracking performance was shown to be substantially improved 
by use of down-linking of velocity information (30/c for position and 80% 
for velocity). However, wind effects were not investigated. 

Recommendations 

It is recommended that the models of the Mode S tracking system be 
incorporated into a fast-time simulation of maneuvering aircraft, and 
that the resulting data be used to determine quantitatively the improve- 
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merits available for trajectory reconstruction. Also, it is recommended 
t^at actual flight data be taken for an instrumented maneuvering aircraft 
tracked by both Mode S and a precision radar tracker. Such facilities 
exist currently at the FAA Technical Center. This data could then be 
used for trajectory reconstruction and accident analysis research. 

Second, it is recommended that the concept of down-linking of aircraft 
velocity information to enhance the precision tracking available from 
the Mode S system be pursued. This includes investigating the following: 

a) The increased accuracy available for tracking the aircraft 
and generating the aircraft trajectory from the tracking data; 

b) The functional requirements necessary to add airborne measurement 
of airspeed, heading, and sink rate to the down-linked data plus 
the requirements necessary to collect, store, and process this 
data at the Mode S sensor site. 

c) The cost and benedits of expanding Mode S to utilize additional 
down- linked variables. This includes applications such as 
improved ATC tracking, collision avoidance (both with ATAS and 
TCAS) , and cockpit displays of traffic information (CDTI), as 
well as enhanced accident analysis. 
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APPENDIX A 


TRACKING ANAI-YSIS 


The radar tracking accuracy that is available from the Mode S system 
can be broken into the vertical dimension and the horizontal plane. The 
vertical accuracy is dependent upon the processing of the encoded altitude 
which is quantized to the nearest 100 ft. The encoded altitude errors are 
further described in Chapter IV. The filtering of this dimension is 
described in Refs. 29 and 30, and is not repeated here. 

The subject of this appendix is filtering of the data describing 
horizontal position and velocity of the tracked aircraft. It is assumed 
that the basic a-6 filter is used for this purpose, and that the horizontal 
position is described by x-y Cartesian coordinates. This implies that 
range and azimuth measurements of the target position are first transformed 
to equivalent Cartesian coordinate measurements. Then, an a-6 filter is 
used to process the data in each orthogonal axis. 

We are interested in several points in conducting this analysis: 

a) The relationship of measurement geometry (between the tracked 
aircraft and the radar antenna) on the resultant estimation 
accuracy. 

b) The effect of measurement error (Mode S and ATCRBS range and 
azimuth errors) on the estimation accuracy. 

c) The effect of aircraft turn rate on estimation accuracy. 

d) The improvement in estimation accuracy which would result 
from down-linking of aircraft velocity as measured on-board. 

In this preliminary analysis we assume that measurement errors can 
be described as zero mean white noise with constant variances. Acceleration 
is treated as an unknown deterministic constant. A more elaborate evalua- 
tion of error effects using the models presented in Chapter V and a more 
complex filter or smoothing technique requires computer simulation at a 
future time. 
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Performance Comparison of Two Horizontal Filter Configurations 


Performance of two filter configurations is now discussed by using 
the current ATCRBS and projected Mode S SSR sensor data error values. These 
two configurations are (a) the non-aided (i,e.> the surveillance data only) 
mode and (b) the aided mode with down- linked aircraft dynamic variables. 

The latter is motivated by the potential future ability of the aircraft to 
communicate with the ground sensor through the discrete addressing channel 
of Mode S. 

The horizontal tracking performance configuration is made ir. the 
following by using generalized dimensions. The estimation errors are 
discussed from a filtering rather than smoothing point of view. There are 
two reasons for this: (a) This is only a preliminary analysis, and (b) 

the error assessments are easier in the filtering context. However, the 
general conclusions from this view point extend to a smoothing approach, 
because the filter algorithms would be the basis of the smoothing process. 


Aircraft Motion and Measurement The aircraft kinematic equations of 
motion are given by the second order equation, 




"l 

L~ 



— 

.0 

1. 

X + 

n 

. A . 


OIF. mok QbrniTY (A.i) 


Here, = state vector of position and velocity at time nL , 

” X 

n 

X 

^ n _ 

acceleration component due to a turn maneuver; and 

surveillance data sampling period (nominally, 4-5 sec 
for the Mode S sensor in the terminal area) . 


a = 
n 

L 


Note here, that we are presenting the filtering problem in only one 

dimension, where in fact, the horizontal plane is defined by two orthogonal 

dimensions. This simplification is strictly for convenience. The results 

can be easily expanded from the two-variable state (x , x ) to the four- 

n n 

• « 

variable state (x ,v ,x ,v ) . 

n ■ n n ' n 
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The surveillance (measurement) data are assumed to be given by 


X 


mn 


= X + 5 
n ^n 


(A.2) 


This ass'jmes that the radar measurements of range and azimuth are resolved 
into appropriate Cartesian coordinates. 


Also, in Eq. (A.2), 5^ is assumed to be an independent white noise 
sequence. For any radar based system, the measurement error consists of 
two parts - range error and azimuth error. Therefore, the variance of the 
measurement noise would be of the form 


o 




2 



+ r^o. 




where is the range error standard deviation, r is the range, and Og is 
the bearing error standard deviation. For the Mode S sensor, the range 
error is relatively negligible. Therefore, 


is a good representation. The model Eqs. (A. 1) through (A. 3) contain 
sufficient characteristics and information to analyze the problem at hand. 
Again, this assumes a one-dimensional problem that can be easily expanded 
to the two horizontal dimensions. 


Non-aided Filter Configuratlofi The non- aided track-while-scan filter 
algorithm is given by the following so-called (two- variable) alpha-beta 
tracker equations: 



n+1 


1 A 
0 1 



a 

S/A 


mn+1 


(Prediction) , 


(A. 4) 


“i\p} 


(Mearurement Update) . 


It is noted that the above algorithm is based upon assuming from Eq. (A. 1) 
that acceleration a^ equals zero. When the acceleration is constant, then a 
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three-variable filter can be used (the so-called alpha-beta-gaimna tracker). 
However, there are two reasons for not using sucli a three-variable filter. 
First, the tracked aircraft longitudinal acceleration (typically, O.lg 
max) would be small compared to the acceleration due to a turn maneuver 
which could be up to 0.5g for an airplane. Secondly, the acceleration 
duo to a turn is not constant along any one dimension. 

The planned Mode S ARTS III tracker algorithm incorporates a three- 
state variable filter. However, the third state "estimated acceleration" 
is used for an ad hoc turn delt^ction logic to compensate for a detected 
turn. The performance of this so-called threshold ai:’> tracker (TABG) is 
analyzed in Ref. 31. Here, we are primarily interested in evaluating 
effects of variable measurement error magnitude and down-linking of 
addi tionai information . 


Now given the filter configuration given by Eq. (A. A), the usual Kalman 
filter technique can be invoked to obtain the optimum gains if the model 
acceleration input is an independent white noise sequence of known intensity. 
However, the primary acceleration is assumed to be the result of a turn 
maneuver rather than a high frequency disturbance due to gusts or pilot 
induced motion* Therefore, the acceleration input must be treated as deter- 
ministic but unkno\^m. In such a case, the filter gains are tuned by 
optimizing a combination of deterministic and stochastic errors assuming 
worst case magnitudes. 


The gain optimization process can be facilitated further, if the gains 
a and S are expressed by a single parameter. t’or this reason, we confine 
this analysis to the critically damped gain configuration. With this assump- 
tion, the gains are given by 


a = 


B - 



(A. 5) 


Here, the parameter is related to the filter bandwidth (= 
the expression 
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Y = exp(~Wj^A) = exp(-A/tj^) . 


ommv m::? 


(A.6) 


With this form for the gains, the estimation error equation is given by 


n+1 


1 

<3 

CM 

< 

CM 

\ 

X + 

2,2 
Y t 

2 

a + 

^ 0 
1~Y 

Y(2-y) 

n 


n 

(1-y)^ 

A 


Here, the estimation error is defined as 


°n+l 


(A. 7) 


n 


X - X 

n n 


(A. 8) 


Equation (A. comes from combining Eqs. (A.1)-(A.6). 


If we assume that the acceleration a is a constant, the above 
equation has a mean steady-state error given by 




a-y)^ 

1h-3y a 

L ■ 1-Y 2 


(A. 9) 


If the measurement noise variance is stationary for a reasonable time 
duration, Eq. (A. 7) can be used to obtain the stationary error covariance 
matrix : 


Cov (x^) = E 


2 2 


|(x x^)!* , 


(A. 10) 


2. 


(1~y) (1+4y+5y ) (1~y) (5+4y+y )/4 


(l-Hy)'’ 1_(1 -y) (5+4y+Y )/4 


2(1-y)^/A^ 


The root-mean-square errors are given by combining the mean and 
variance terms from Eqs. (A. 9) and (A. 10). This gives the expressions. 
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Y^a'^ 1 (1-y) (1+4y + 5y^) 2 2 

4 ® 3 ^ “b 

(1-y)^ (1+y)^ 


(A.ll) 


(1+3y)^ a! ■ r^a 2 

(1-y)^ ^ (n-y)^ A ^ 


Each of these expressions depend on the five parameters: 
a^ sampling rate, L ; 

b) acceleration magnitude, a ; 

c) measurement error magnitude, a ; 

p 

d) range , r ; and 

e) the filter bandwidth, y ~ exp(-oj,A). 

b 

The first four can be used to obtain the optimum bandwidth. 


As an example of optimizing the bandwidth, the following Mode S 
parameter values represent a worst possible case: 

L = 4 sec , 


= 0.04^ 


= 60 nmi. 


= 255 ft. 


Three values of bandwidth - 0.115, 0.175, and 0.267, taken from Ref. 32, 

were used in Eqs (A. 9) , (A. 10), and (A.ll). Table A.l shows the magnitudes 

of various errors under the assumed conditions. In Table A.l, it is seen 

that ^o liinimize rms position error, the bandwidth u), should be at the 

b 

0.267 value. To minimize velocity error, the higher bandwidth is also 
better. 

The rms errors given by Eq . (A.ll) are very general and applicable 
to many surveillance tracker problems which would fit the basic assumptions 
given earlier. Two applications are considered below. One is the 
sensitivity of errors with respect to range. The other is the effect of 
the bearing error. The currently operational Air Traffic Control Radar 
Beacon System has a bearing error of approximately 0.25 deg (Ic). This 
is considerably larger than the proposed Mode S sensor, 
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Table A.l. Non-Aided Filter Configuration Errors (Mode S) 


bandwidth 

•^b 

rad/sec 

Position Errors 

(ft) 

Velocity 

Errors 

(i 

ft/sec) 

(Eq.A.9) 

a 

(Eq.A. lO: 

X 

rms 

(Eq.A.ll) 

E(x^) 

(Eq.A.9) 

a* 

x^ 

(Eq.A. 10) 

• 

X 

rms 

(Eq.A. 11) 

0.115 

750 

180 

771 

251 

10 


251 

0.175 

249 

209 

325 

158 

25 


160 

C.267 

1 

i 

70 

235 

245 

99 

32 


104 


Table A. 2 shows the nns error at several range values for Mode S and 

ATCRBS system characteristics. It was assumed that the surveillance period 

for both was 4 sec. The rms errors are given at the filter bandwidth of 

0.267 rad/sec. Therefore, it should be noted that the rms errors based on 

optimized filter gains would be smaller than the table values. Also note 

that the speed errors (x ) are dominant by the acceleration due to turn 
- ^ rms ■' 

term (x = 99 ft/sec), 
a 


Table A. 2. Geometric (Range Magnification) and ’ 

Sensor Bearing Error Effects 



(■^b ~ 

0.267 rad/sec) 


— 

Range 

Mode S (a_ = 0.04) 



ATCRBS (a. =0.25) 

p 

(nmi) 

S (ft) 

rms 

X (ft/sec) 

rms 

rms 

X (ft/sec) 

rms 

15 

91 

99. 

374 

111 

30 

137 

100. 

738 

141 

45 

190 

102. 

1104 

179 

60 

245 

104. 

1470 

i 

223 

1 _■ - ■■ - - - . -- — _| 


Aided Filter Configuration For the aided filter configuration with 
down-linked aircraft variables, it is assumed that a measure of ground 
velocity is available from airborne measured true airspeed and heading 
information. In our previous study [6] , an approach was taken to obtain 
a turn rate estimate from the heading information using a two state variable 
filter. Then, the resulting turn rate estimate was used to derive position 
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and vDlocicy estimates. This approach was taken for two reasons: One 

was that the system transition matrix contains functions of turn rate, 
and the other was that the horizontal components have strong cross-coupling. 
(The latter implies that the dynamics due to turn maneuver cannot accurately 
be described by a single dimension.) Never-the-less , tne above approach is 
equivalent to obtaining an estimate of the average acceleration over the 
last sampling period using the down-linked velocity information. Thus, the 
observed acceleration a is assumed to be of the form 


Here , 


a 

n 

a 


n 


+ n. 


(A. 12) 


observed average turn acceleration, 
true acceleration, 

acceleration observation error (again assumed to be 
an independent white noise sequence) . 


The error r may be thought of as consisting of two components. One, o^, 
is due to the measured true airspeed error magnitude, and is due to 
the measured heading error magnitude. In this simple analysis, we are 
ignoring the effects of winds and the fact that measured true airspeed is 
not groundspeed. Thus, 


c 


2 


k. a 
1 V 


1 2 2 

■"a 


> 


where and are attenuation factors. In this preliminary study , it is 
assumed that 

2 

a = 1.5 kt/sec = 2.5 ft/sec 

n 


With this additional downlinked data, the aided tracker configuration 
is given by the following equations. 



’1 L~ 



"■a^/2" 

X 

np 

.0 1_ 
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_ A _ 



[“1 
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^n+1 

X + 

np 

r" 

u> 

t> 

1 _ 


mn+1 


5 (Prediction) , 

(A. 13) 

[J. Olx^ i (Me asurement Update) . 


Here, the same notation (a and 6) is used for the filter gains. 




With the same critically damped gain configuration of Eq. (A. 5), the 
es imation error equation is given by 


'n+1 




v(2-.) 




2,2 

I_A 

2 





1 

CM 

1 

n + 


n 

(l-Y)^ 


A 


’n+1 • 


(A. 14) 

Equation (A. 14) is similar to Eq. (A. 7) except that Eq. (A. 14) contains 

the white noise term n rather than the deterministic term a_. 

n n 

Assume that the errors are zero mean, independent, white noise sequences; 
the mean error is zero, i.e.. 


E {x } = 0 . 

n 


(A. 15) 


Thus, in this configuration, there is no dynamic delay error due to 

unmodeled turn acceleration. Under the assumption of stationary statistics, 

the steady-state error covariance (P = cov(x)) satisfies a linear equation 

s s 


T T 2 T 2 

P - F P F a ^ 

ss cssc lln 22? 


(A. 16) 


Here, F^, , and G 2 are the system matrices from Eq . (A. 14). Again, the 

covariance matrix depends on five parameters. In order to obtain a good 
value for filter bandwidth by minimizing the covariances, Eq , (A. 16) was 

solved numerically for several values of the parameter, . Table A. 3 shows 

the results of this computation. The magnitudes of , Yj and A were 

as before. For the cases chosen, the lower error values for range of 60 nmi. 

a = 170 ft. 

X 


a* = 

X 


19 ft/sec. 


at the filter bandwidth of 0.09 rad/sec. For 15 nmi. range, these reduce 
to 50 ft and 12 ft/sec with a bandwidth of 0.25 rad/sec. 
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Table A. 3. Aided Filter Configuration Errors at 60 nmi. 


Bandwidth 
u, (rad/sec) 

D 

Position Error 
0 (ft) 

X 

Velocity Error 
O' (ft/sec) 

X 

0.05 

220 

25 

0.09 

170 

19 

0.11 

176 

19 

0.13 

185 

19 

0.17 

201 

21 


Conclusions By comparing the values in Tables A.l and A* 2, it is seen 
that the large velocity errors preclude using estimated velocity for any 
precision applications. For both Mode S and ATCRBS, the range magnification 
effect is larger for the position errors than the velocity errors. This 
means that the velocity errors are primarily due to the dynamic delay effect 
of the unknown acceleration. The bearing error of ATCRBS (0.25 deg) has 
a significant effect on the estimation errors (0.06 - 0.24 nmi position error) 
when compared to the Mode S values based on 0.04 deg. It should be noted, 
however, that the gains are not optimized. Another notable point is that 
the deterministic and stochastic errors are weighted equally in Eq . (A. 11), 

i.e., this assumes that the turning and the straight flight portions are 
equally important. In reality, most of a flight trajectory consists of 
straight segments (connected by short circular arc segments). Therefore, 
gains may be chosen so that maneuver tracking performance is degraded in 
order to yield a better performance in straight segment tracking. 

By comparing the values in Tables A.l and A. 3, it is seen that the 
(down— linked variable) aided filter^ s bandwidth is much narrower (0.09 vs. 
0.267 rad/sec). Furthermore, the position error is improved by about 30?:, 
and the velocity error is improved by about 80%. This preliminary analysis 
shows that the aided filter configuration would improve the state estimates 
substantiallv when compared to the non-aided filter configuration. More 
detailed analysis is required to assess the effects of winds. 
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Ont can conjecture that comparable conclusions would apply to the 
smoothing process. This is because the smoothed estimate is a linear 
coroi .ation of forward and backward filter estimates. Thus, the filter 
ptrfjrmance essentially dictates the smoother performance. Therefore, the 
limiting factor would be the basic accuracy achievable by a particular 
filter configuration. The preliminary analysis shows that the non-aided 
filter configuration has a basic limitation when the tracked aircraft is 
turning . 

Smoother performance was not analyzed because its error characteriza- 
tion is complex. However, it would be interesting to perform this analysis 
using realistic simulated data at a future date. 


Tracking Errors for Orbiting Targets 


The previous section examined the effects of measuremei.t errors and 
the advantages of down-linking of velocity information to the ground 
tracker. Errors were characterized as white noise or as a constant unknown 
acceleration due to turning. In this section, we examine the expected 
tracking errors using a different method. Here, the tracked target is 
assumed to be flying in circles around the tracker (i.e., it is orbiting.) 
This causes sinusoidal oscillations in values of the Cartesian position 
coordinates. The oscillations enable us to examine the filter* s perfor- 
mance by using frequency analysis techniques. 

These results are a summary of the analysis found in Ref. 32. In 
Ref. 32, it was assumed that a-6 filters were used in an on-board tracking 
system to process range/bearing data for a collision avoidance system. 

The target aircraft was assumed to be flying in orbits around the CAS- 
equipped aircraft. However, these results are directly applicable to 
ground-based radar tracking. Figure 1 describes t)»e simple two-dimensional 
tracking system making use of identical ct-6 filters in each of the two 
Cartesian coordinates. Equations (A. 4) are again applicable for describing 
the filters. 
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The method of z-transforms is used to obtain the discrete transfer 

function between the smoothed estimates of x and x„, the predicted 

n n 

estimate x , and the measurement x . These transfer functions are 
np mn 


X 


z(az + (S-g)) 


z" - (2-a-S) 


-o-8)z + 1-0 Poor 


QUALi 


S3 

ty 


(A.17) 


H’ 


£z ~ 8z 


(z - (2-a-0)z + l-a)A 


(A. 18) 


np 


(g+ g) z -g 


z - C2-g-6)z + 1-g 


(A. 19) 


Here, z is a complex variable. The frequency response is obtained by 
letting z = exp(ju)A) with ui in rad/sec. 


The ability of the g-8 filter to smooth the uncorrelated input 
measurement noise is characterized by the variance reduction ratio. 
This is defined as 


noise variance (power) of response signal 
noise variance (power) of input signal 


(A. 20) 


In Ref. 32, this is found to be 


K, = |h(z)(z-Zt) H(z ^)/zl +iH(z)(z-z.,) H(z ^)/zl .(A. 21) 

( 1 I z-z^ ( z , ^”^2 

Here, z^ and Z 2 are the poles of the g-6 filter. 


The poles of Eq. (A. 17) -(A. 19) are found to 

^1,2 " 1 ±y /(a+E)^ - 4B 

The poles are substituted into Eq. (k.22) to yield 


(A. 22) 
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(A. 23a) 
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(A.23b) 


(A. 23c) 


Equations (A. 23) relate the a-g filter output variance to the input 
measurement variance. 


a 


Figure A. 2 shows the target aircraft moving in a constant radius, 
constan" speed orbit. Even in the absence of measurement errors, the 
a-£ filter produces tracking errors as discussed before. These are 
derived to be 


E {x i 
a 

E {x } 
a 


R (H- - 1) , 

ox 

n 

- J") • 

n 

R (H - 1) . 

o X 

P 


(A. 24a) 
(A.24b) 
(A.24c) 


where R^ is the orbit radius and u- is the target angular rotation rate. 
Also, from Fig. A. 2, 


w = a/v , (A. 25) 

a = v^/R^ . (A. 26) 


where v is the aircraft speed and a is the associated acceleration. Note 
that the tracking errors confuted from Eqs. (A. 24) are complex. They 
indicate magnitude and direction as depicted in Fig. A. 2. 


The magnitudes of the tracking errors given by Eqs. (A. 24) are 


approximated 

by 




!e 

{x }; 
a ' 


(l - a)aA“/g , 

(A. 27a) 
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(A. 27b) 
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(A. 27c) 
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Figure A. 2. Geometric Interpretation of Orbiting Target 

Tracking Errors [321. 



The validity of these approximations is illustrated in Fig, A. 3) for 
acceleration set at 0.5 g, ^ of 1 sec, and (a,£) of (0.4, 0.1), respectively. 

Values of a and 6 which minimize noise output variance tend to 
increase tracking errors for a tuiming target, and vice versa. One method 
of choosing a and B values is to minimize the root-simi-square (rss) values 
of the total error from Eqs. (A. 24) and (A. 27). In Ref. 32, these values 
were numerically solved for 0.5g acceleration, o of 825 ft and L of 
1 sec and 4 sec. The optimum values to minimize the position estimate 


error were found to be: 



4, (sec) 

a 

6 

1 

.25 

.066 

4 

.60 

.431 


Reference 32 lists -he Fortran code to find these values. Figure A. 4 shows 

the composite errors for various values of a and 6 with different target 

speeds. Note that these values ccii^are with those shown in Table A.l if 

the larger value of a is taken into account. As can be seen, there is 

xm 

no unique pair of (a, 6) values that minimizes all three error terms. Other 
values of acceleration, measurement error, and sampling interval will 
result in different optimum values of a and B. 
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